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Abstract. The course of 5 lectures given at the seminar "Integrable Systems: from 
Classical to Quantum" (Universite de Montreal, Jul 26 - Aug 6, 1999) contains a de- 
tailed comment on the recently discovered (Gaudin-Pasquier, 1992) connection between 
Backlund transformations in the theory of classical integrable systems on one hand, 
and Baxter's Q-operator for quantum integrable systems, on the other hand. We re- 
, strict our attention to the systems with finite number of degrees of freedom. Our main 

Q^ I illustrative example is the periodic Toda lattice. We present a general construction 

O ■ of Q-operator for models governed by the 5'L(2)-invariant i?-matrix and apply it to 

our example. We discuss also applications of BT and Q-operators to the separation of 
O ■ variables and theory of special functions. 
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1 Introduction 



1.1 A bit of history 

The Backlund and Darboux transformations appeared in XIX century in the study 
of the problems of differential geometry. With the advent of the Inverse Scattering 
Method in 1960s, their relevance to the integrable nonlinear evolution equations was 
quickly recognized, and the amount of literature accumulated since then is enormous. 
See, for example, the monograph |]2^. Especially important for these lectures is the 
Hamiltonian interpretation of BT discovered by Flaschka and Mclaughlin \T2\. 

The Q-operator belongs to the realm of quantum integrability, and, compared to 
BT, is a relatively new invention. The Q-operator was introduced first by Rodney Bax- 
ter in his seminal study 0, see also [§], of the integrable quantum XYZ spin chain as an 
ingenious device which allowed to determine the spectrum of the model — the problem 
which was intractable by other known methods, such as Bethe ansatz. The Q-operator 
is actually a one-parametric family of operators Qx commuting with the Hamiltonians 
of the integrable system. Its main characteristic property is that its eigenvalues satisfy 
certain finite-difference, or, depending on the integrable model, differential equation 
with respect to the parameter A, known nowadays as Baxter equation. The Baxter 
equation, together with appropriate boundary conditions, provides a one-dimensional 



multiparameter spectral problem which allows to determine the spectrum of the com- 
muting Hamiltonians of the model in question. Thus an originally multidimensional 
spectral problem is reduced to a one-dimensional one — phenonomenon similar to the 
separation of variables (SoV). The coincidence is not an acccident. Indeed, as shown in 
19, ^ , for classical Hamiltonian systems there exists an intimate relation between SoV 



and Backlund transformation (BT), the latter being the classical analog of Q-operator 



19, 25 



For a long time the XYZ model remained the only model for which a Q-operator 



was known. In 1992 Pasquier and Gaudin p5[ have constructed a Q-operator for the 



quantum periodic Toda lattice using a somewhat different approach from Baxter's one. 
They have described Q-operator explicitely, as an integral operator, and have found an 
important relation between the Q-operator and the Backlund transformation from the 
classical Toda chain. Namely, the Backlund transformation, as a canonical transforma- 
tion, coincides with the classical limit h ^ oi the automorphism O i-^ Q\OQ^^ of the 
associative algebra of quantum observables generated by Q\. The generating function 
Fx{y\x) of the Backlund transformation is obtained from the semiclssical asymptotics 
Q\{y\x) ~ ey.\){ifi~^ F\{%)\x)) of the kernel Q\{y\x) of Q\ considered as an integral 
operator in the coordinate representation. 

Later on, in P] Bazhanov, Lukyanov and Zamolodchikov gave a boost to the original 
Baxter's idea of constructing Q\ as the trace of the monodromy matrix constructed of 
Lax operators corresponding to a special representation of the relevant quantum group 
in the auxiliary space. Setting the problem in the context of representation theory 
for quantum groups, they have taken as such representation of slq{2) the so-called q- 
oscillator representation and have managed to construct a pair of Q-operators for the 
massless sine-Gordon quantum field theory in a periodic box. It seems that the same 
g-oscillator representation allows to construct a Q-operator for any integrable model 
governed by the quantum group slq{2), see |Q. 

In the paper |2^ devoted to construction of a Q-operator for the so-called DST 
(dimer self-trapping) model (a degenerate case of XXX magnetic chain) the combi- 
nation of the approaches due to Baxter (Qa as a trace of monodromy) and Pasquier- 
Gaudin {Qx as an integral operator) allowed to describe the structure of the Q-operator 
in the greatest detail. Besides an explicit expression for the kernel of the integral op- 
erator Q\ in several equivalent forms one can calculate explicitely the matrix elements 
of Q\ in the natural monomial basis. 

In the last years a considerable progress is achieved in the understanding of the 
Hamiltonian properties of the BT for the classical integrable systems which parallel 
those of the Q-operator in the quantum case. It is worth noticing that the classi- 
cal counterparts of some of the properties of the quantum Q-operator were unknown 
before. As an example one can mention the so-called spectrality property of BT dis- 
covered in [|T9] which corresponds to Baxter's finite-different equation for the quantum 
case. Baxter's construction of Qx as a trace of monodromy matrix has led to a new 
construction of BT from symplectic leaves of the quadratic r-matrix Poisson bracket 

A special topic actively studied in the last years is the relation of BT and SoV which 
turns out to be twofold. On one hand SoV can be obtaind from a composition of BT or 



Q, see [O. On the other hand, a BT can, in turn, be obtained as the transformation 



intertwining a pair of SoV [^. In this case, the quantum interpretation is not found 
yet. 

The growing interest in studying various properties of Q-operator for a variety of 



quantum integrable systems is indicated by a surge of recent pubhcations, see |T^, ^, 

m- 

1.2 Plan of lectures 

In these lectures I will concentrate on the parallels between the Backlund transfor- 
mation for the classical Hamiltonian sytems one one hand and the Q-operator for the 
quantum integrable systems, on the other hand. As it frequently happens, when two 
theories merge after having been developed for a considerable time independently, the 
resulting cross-fertilization is quite useful for both. The most recent example is the 
classical r-matrix and Lie-Poisson groups |2^, |3^ whose invention was inspired by the 
quantum theory. 

We shall restrict our attention to the systems with a finite number of degrees of 
freedom (pure quantum mechanics, no field theory) and put special stress on Hamilto- 
nian mechanics which is essential for quantization. All the new notions and techniques 
will be introduced on the example of the periodic Toda lattice and accompanied with 
a short discussion of possible generalizations. 

The lectures can be considered as an extended commentary to the paper by Pasquier 
and Gaudin p^ accompanied by the original results obtained by V. Kuznetsov and 
myself pi, EUl 



2 Classical periodic Toda lattice 

2.1 Description of the model 

The periodic Toda lattice ||38|, |23| is a system of n degrees of freedom described in terms 



of canonical coordinates x = (xi, . . . , x„) and momenta X = (Xi, . . . , X„) having the 
standard Poisson brackets 

{Xj, Xk} = {xj, Xk} = 0, {Xj, Xk} = Sjk. (1) 

In what follows we always denote canonical coordinates with small letters, e.g. x, 
y, s, t, (f and the corresponding canonical momenta with the respective capital letters: 
X, Y, S, T, $. Such a convention helps to deal with several sets of canonical variables. 

The physical Hamiltonian H of the Toda lattice 

^ = EQ^I + e'^"^"'^) (2) 

describes the system of n one-dimensional non-relativistic particles of equal mass in- 
teracting via exponential potential between the nearest neighbors. In formulas like @ 
we always assume the periodicity convention: j + n = j. The Hamiltonian H is thus 

4 



invariant with respect to the translation j ^^ j + 1, hence the name 'periodic Toda 
lattice'. 

The equations of motion / = {H, /} corresponding to the Hamiltonian (H) are 

Xj = Xj, (3a) 

Xj = -e^i-^i-i + e^^+i-^^ (3b) 

2.2 Integrability 

The periodic Toda lattice is an example of a completely integrable Hamiltonian system 
in the Liouville- Arnold sense f^. It means that the Hamiltonian H (0) is an element 
of a ring generated by n independent Hamiltonians Hi,. . . ,Hn which commute 

{H,,,H,,} = (4) 

with respect to the Poisson bracket (|l|). As a consequence, the Hamiltonian flow 
/ = {H, /} leaves the Hamiltonians Hj invariant Hj = and, therefore, leaves invariant 
the level manifolds Ph obtained by fixing the values of the Hamiltonians Hj = hj. 

The fundamental result in the theory of Hamiltonian integrable systems is the Liou- 
ville theorem [Q which claims that the level manifolds Ph, if compact, are diffeomorphic 
to n-dimensional tori T". Moreover, there exist such canonical action-angle variables 
$, (f that the action variables ^j are functions of the Hamiltonians Hi,. . . ,Hn, and the 
Hamiltonian flows linearize in the angle variables: {Hj,ipk} = ^jk{^)- 

The easiest way to demonstrate the integrability of the periodic Toda lattice is to 
make use of the Inverse Scattering method (or, the Isopectral Deformation method) 
in its Hamiltonian version, see for example [|TI| and Prof. J. Hamad's lectures in this 
volume. Within the ISM framework, the commuting Hamiltonians are obtained from 
the spectral invariants of the Lax matrix L{u; X, x) which is a square matrix of order 
N (generally speaking, different from the number of degrees of freedom n) depending 
on a complex parameter u called spectral parameter and whose matrix elements are 
functions on the phase space. 

The spectral invariants tk{u) of L{u) defined as the coefficients of the characteristic 
polynomial 

N 

W{v, u) = det(t; - L{u)) = w^ + ^(-l)'=i;^~''tfc(u) (5) 

fc=i 

are elementary symmetric polynomials of the eigenvalues of L[u) , or, in terms of matrix 

elements of L{u), sums of principal minors of order k (determinants of submatrices of 

L[u) of order k whose diagonal is contained in the diagonal of L{u)). For example, 

ti{u) = trL{u), t]\f{u) = {—l)^detL{u). The commuting Hamiltonians Hj are usually 

obtained as coefficients of tk{u) when tk{u) are polynomials in u, or coefficients of 

expansions of tk{u) in some other bases of functions of u (e.g. trigonometric or elliptic 

ones). 

Proving the commutativity (^) of the Hamiltonians Hj is equivalent thus to proving 

the commutativity 

{'tkiiui),tk2iu2)} = (6) 



of the spectral invariants tk{u). The following theorem due to Babelon and Viallet 
provides a technical mean to do it. As proven in |0, the the commutativity (P) of the 
spectral invariants of L{u) is equivalent to the existence of a so-called r-matrix repre- 
senation for the Poisson brackets {Laibi(Mi), //a2 62('"2)} between the matrix elements 
Lab{u). To write down the representation in a compact form, we introduce the tensor 
product notation 

L = L O 1, L = 1 O L, (7) 

1 2 

where 1 is the unit matrix of order A^. Respectively, {L (ui) , L {U2)} is the matrix of 
order A^^ x A^^ of all Poisson brackets between the matrix elements of L{ui) an L{u2)- 
The theorem of Babelon and Viallet claims that the commutativity @ is equivalent 
to the existence of two matrices, ri2 and r2i, of order A^^ x A^^ such that the equality 

{L{ui),L{u2)} = [vi2,L{ui)] - [t2i,lM] (8) 

holds for any mi, M2. Note, that, X\2 and X2\ depend on u\ and ^2, and, generally 
speaking, contain the dynamical variables X, x. Actually, one can always choose the r- 
matrices in such a way that t2i(Mi, M2) = "^12^12(^25 ^i)'Pi2 where V\2 is the permutation 
matrix: ^123; ®y = y®x. 

Speaking again about the periodic Toda chain, in order to construct the commuta- 
tive Hamiltonians, we have to produce a Lax matrix and the corresponding r-matrices. 
There are at least two possible Lax matrices for the periodic Toda chain, one of order 
2x2, another one of order n x n, see |Tl|, |3^, |23|, |l|. In this subsection we shall work 



with the 2x2 matrix. 

The 2x2 Lax matrix (or, monodromy matrix |lT|) L{u;X,x) is defined as the 
product of local Lax matrices ij{u) depending on the variables Xj and Xj only: 

L{u)=£n{u)...e2{u)ei{u), (9) 

i,{u) ^ i,{u; X,, X,) = ( ""t^^ -f \ . (10) 







The characteristic polynomial of L{u) is quadratic in v having thus two spectral 
invariants: ti{u) and t2{u). However, t2{u) = det L{u) = 1 by virtue of det i{u) = 1 
which leaves t{u) = ti{u) = ti L{u) as the only nontrivial spectral invariant: 

W{u, v) = det(i; - L{u)) = v"^ - t{u)v + 1. (11) 

The Hamiltonians Hj are obtained then from the expansion of t{u): 

t{u) = u" + ifiu"^^ + ... + Hn. (12) 

In particular. Hi = Xi + . . . + X„ is the total momentum. It is easy to see that the 
physical Hamiltonian (^ is given by the formula H = ^Hf — H2 and thus belongs to 
the polynomial ring generated by Hi,. . . ,Hn- 

To prove the commutativity (^) of the Hamiltonians Hj, or, equivalently, the com- 
mutativity 

{t{ui),tiu2)} = (13) 

6 



of their generating function t{u), it is sufficient to find the corresponding r-matrices. 



Actually, what we are able to prove is a much more special representation |jTT| for the 
left-hand-side of (E): 



where 



{Liu),Liv)} = [n2iu-v),Liu)L{v)] (14) 

r^iui - U2) = ^-^— (15) 

Ml - U2 

is the S'L(2)-invariant solution to the classical Yang-Baxter equation 

[ri2{u),ris{u + v)] + [ri2iu),r23{v)] + [ris^u + v),r23{v)] = 0. (16) 

One can easily transform the formula (0) to the Babelon-Viallet form by setting 
m(i 

1 22 

^i2 = -{ri2{ui -U2)L{u2) +L{u2)ri2{ui -M2)), (17a) 

1 1 1 

t2i = -2(^12(^1 -U2)L{ui) + L{ui)ri2{ui -U2)). (17b) 

To prove (|1^ we shall make use of the so-called comultiplication property of the 
quadratic Poisson bracket. It is a simple exercise to verify that if two Lax matrices 
Li{u) and L2{u) defined, respectively, on different phase spaces Pi and P2 satisfy each 
the identity (|I4D then their matrix product L{u) = Li{u)L2{u) defined on the direct 
product Pi X P2 satisfies the same identity. The proof uses nothing but the identity 



(14) for Li and L2, and the identity 

{ii(ni),l2(n2)} = 0. (18) 

It is sufficient thus to verify the identity ([T^) for the local Lax matrices ^j{u) given 
by (^) which is a matter of a direct calculation. 

Strictly speaking, to establish the integrability, besides the proof prove of the com- 
mutativity of Hamiltonians Hj we need to prove their independence. For a proof, see 
p8| . It is also possible to verify that, modulo center-of-mass motion, which is easily 
separated, the level manifolds P^ are compact and thus, by virtue of Liouville theorem 
are isomorphic to tori. 

2.3 Quadratic Poisson bracket 

Before starting the discussion of Backlund transformations we need to learn some more 
facts about the r-matrix quadratic Poisson bracket ([I^) and the class of integrable 
models it generates. 

Let us suppose that 2x2 matrix L{u) is a polynomial in u of degree n 

L{u; X, x) = L(")m" + L("-i)m"~i + . . . + L^°^ (19) 



and regard the equality (JT^) with the r-matrix given by (|T3p as introducing a Poisson 
bracket on the 4(n + 1) variables L|fj, , j = 0, . . . ,n, a,b = 1,2. For the sake of simplicity 
we think of L^^jj as complex variables and do not consider here the question of choosing 
an appropriate *-conjugation. 

It is easy to see, that, despite the denominator [ui — U2) present in the r-matrix 
([T5|) , the right-hand-side of (14) is polynomial both in ui and U2 because of the identity 



["P, L ® L] =0 which nullifies the numerator for ui = U2- According to a theorem by 
Sophus Lie |^, for any Poisson bracket there exist local coordinates (X, x, c) such that 
X and X are canonical (|l]), and c are central, that is 

{cj,Ck} = {cj,Xk} = {cj,Xk} = 0. 

To obtain a symplectic manifold which can serve as a phase space for a mechanical 
system, one needs thus to restrict the Poisson bracket onto a level manifold of its 
central (or Casimir) functions. 

In case of the bracket (0) the Casimir functions can be found easily. First, the 
leading coefficient L'-"-' provides 4 casimirs. More casimirs are given by the coefficients 
of the determinant detL(M). Being, generally speaking, a polynomial of degree 2n, the 
determinant has {2n + 1) coefficient but its leading coefficient coincides with det L^"'^ 
which gives us only 2n new casimirs. In total, we have {2n + 4) casimirs which cor- 
responds to the level manifolds of dimension 4(n + 1) — {2n + 4) = 2n. To show that 
there are no more casimirs, it is sufficient to construct an example of 2n- dimensional 
symplectic leaf of the bracket (0). 

The tool for constructing such examples is the comultiplication property of the 
bracket ([l^) mentioned in the subsection |2.2| . It allows to build multidimensional 
symplectic leaves from simpler blocks. The simplest, 0-dimensional simplectic leaf of 
the bracket (|l^ is given by a constant matrix L{u) = K. The next most natural choice 
is to take a linear polynomial in u with the unit matrix as the leading coefficient: 

e^^^iu) = «]i + §, ^=[s, + iS2 ^'-sT ) • ^^°^ 

Substituting ( pOD for L{u) into (p^ we obtain for Sa the Poissonian algebra iso- 
morphic to the Lie algebra s/2: 

3 

{Sa, S13} = — i 2^ ^a/S-yS^, (21) 

7=1 

EaiBj being the standard antisymmetric tensor. The Poisson bracket (^) has the 
Casimir function C = Sf + S2 + S^, and its generic symplectic leaves C = const 7^ 
are 2- dimensional spheres. 

Noting that, due to the fact that the r-matrix (|l^) depends only on the difference 
{ui — U2), the shift of the spectral parameter -u 1-^ m — c is an automorphism of the 
poissonian algebra (pjlf). Therefore, taking a direct product of n copies of the triplets 
Sa restricted to the level surfaces C = p'j, j = 1, . . . ,n we obtain a 2n-dimensional 
symplectic leaf of the bracket (0) given by the product 

L^'^^(m) = Kef'^iu -Cn)... ir^'iu - ci). (22) 

8 



Note that the number of parameters contained in the symplectic leaf is is (2n + 4) 
where 4 comes from the constant matrix K and rest from n casimirs pj and n shifts 
j = 1, . . . ,n. The parameters are easily identified with L^") = K and the zeroes 



c,- 



J' 



of the determinant det L{u) = det K llj{u — Cj — Pj){u — Cj + Pj)- We are thus led to 
the conclusion that the constructed symplectic leaf is in fact the generic leaf for the 
bracket {^). 

The Lax matrix (^) defines an integrable system known as the inhomogeneous 
Heisenberg magnetic chain |[TT|, |2^. All other integrable models associated with the 



Poisson bracket (JT^) and the s/2-invariant r-matrix ([T^) can be obtained from degen- 



erations of the Lax matrix (|2^) . 

To describe some important degenerations of p^ ) let us parametrize the spin com 
ponents Sa in (^) using a pair of canonical variables {X, x) : 

- p —x^X + 2px \ 




'^-^'u) ^ ( '^ -r -- - P - - ^ -/^- (23) 

u — xX + p I 



Multiplying i^^^{u) from the right by the diagonal matrix diag(l, — l/(2p)) and 
performing the shift u ^-^ u + p (note that these are legal operations which do not 
change the Poisson bracket (|T^)) we are capable to take the limit p —>■ 00. The result 



is the Lax matrix for the so-called dimer-self-trapping (DST) model [pC 



«DST/ 




«) = ^ 1 • (24) 



Note that the determinant det i°^'^{u) = —u is linear in u. A further degeneration 
of DST model produces the Toda lattice. To this end, one multiplies £°^^('u) from the 
right by the matrix diag(l, a~^) and, after making substitutions -u 1— > -u — a, x 1— >■ ae^, 
X 1-^ e~^(l + a~^X), obtains in the limit a ^ 00 the unimodular Lax matrix (|ToD for 
the Toda lattice. 

More symplectic leaves can be obtained by applying the automorphism £(-u) h^ 
l{u) = i~^{—u) of the r-matrix Poisson algebra (|1^) to £°^'^(m) and i'^°'^'"{u). Up to a 
scalar factor we have: 

^°"(") ~(x u :lx ) • '™-(") ~ ( e-^ u~-'x ) (25) 

(on i^^^{u) the automorphism acts trivially: i^^^{u) ~ i^^^{u)). 

Conjecture. Any symplectic leaf L{u) of the quadratic r-matrix Poisson bracket 
( my which is polynomial in u can be decomposed (in a non-unique way, of course) into 
a product of a constant matrix K and linear matrix polynomials of the form i^^^{u — c), 
(u-c), F'^^iu - c), F-^^iu - c), F-^^iu - c). 



ODST{ 



In case of a generic symplectic leaf the factorization (|22| ) in terms oi i^^^{u — c) only 
should suffice. The difficult part is to analyze the degenerate cases when the leading 
coefficient L^"'' is a degenerate matrix and/or the degree of detL(M) is less than 2n. 
Hopefully, one of the readers will provide a proof soon. 

One can find more information about the properties of the r-matrix Poisson bracket 



(|TJ) in the papers [^ ^, |2£] as well as in the lectures by Harnad and Reshetikhin in 
the present volume. 



2.4 Backlund transformation and its properties 

In this section we start to study the Backlund transformation for the periodic Toda 
lattice. The Backlund transformation B\ depending on a complex parameter A is 
defined as the mapping from the variables (X, x) to {Y, y) given implicitely by the 
equations 

Xj = e^^-^^ + e^^+i-^^ - A, (26a) 

Yj = e'^-y^ + e^^-^^-i - A. (26b) 



The equations (^61) are algebraic in momenta and exponents of coordinates. Re- 
solving (p9a) with respect to e^^+^: 



e 



Vj+i — p^UY ■ ^ \\ — o'^^J-Vi 



(f^{X.j + \)-e'''^-y^ (27) 

and iterating the equation (|27|) for j = 1,2, ... ,n we finally arrive to a quadratic 
equation for e^^ which implies that the transformation Bx is a two-valued algebraic 
function in terms of X, e^. Fortunately, for all our purposes the simple implicit formulas 
(|26|) are sufficient. 

Another attractive feature of the equations ( pBf ) is their locality: they involve only 
the variables with the indices differing by and 1. Note that even for real A resolving 
the equations (p6| ) can produce complex values of Y and y. To avoid complications, 
we shall not make attempt to study the reality conditions and treat both (X, x) and 
(Y,y) as complex variables, in the spirit of algebraic integrability [0]. 

We start the list of properties of the Backlund transformation with noting its canon- 
icity: the variables {Yj, yj) are canonical. It can be seen from the fact that the equations 
(p6|) can be written down in the form 

X - ^^^ Y- ^^' r28) 

where Fx{y; x) 

n 

Fxiy, x) = E(e"^"'^ - e^^^^"'^ - A(x, - %■)) (29) 

is the generating function [|5[ of the canonical transformation. 
The next property is the invariance of Hamiltonians: 

H,{X,x)=Hj{Y,y), j = l,...,n. (30) 

Though the invariance of physical Hamiltonian H (^ can be proved by a direct 
calculation |]38|, |T^ , to prove the invariance of the whole set of commuting Hamiltonians 
Hj we will need some more effective technique. The easiest way is to make use of 
Inverse Scattering Method explained in section |2]^. The invariance of Hj under Bx 
is equivalent then to the invariance of the spectrum of L{u) which implies that there 
exists an invertible matrix M{u, A) such that 

M{u,X)L{u;Y,y) = L{u;X,x)M{u,X), Vm G C. (31) 

10 



Such a matrix is called Darhoux matrix, and the tranformation of L given by (^TJ) 
is called Darhoux transformation [^ . 



In our case, due to the factorization (y) of L{u) into local Lax matrices ij{u), we can 
be more specific about the structure of Darboux transformation. Setting Mi{u,X) = 
M{u, A) we introduce matrices Mj+i{u, A), j = 1, . . . , n — 1 inductively as 

M,+i(m, A) = ijiu; Xj, x,)M,{u, \)q\u; Yj, y^) 

(note that for j = n, due to the periodicity n + 1 = 1 we recover the equality (|3TD). 
The global transformation (^) takes thus form of the local gauge transformation 

Mj+^{u, \)ij{u; Yj, y,) = ij{u; X„ Xj)Mj{u, A). (32) 

The converse is also true: from (|32|) it follows that the spectrum oi L{u) is preserved. 
To prove the invariance of Hamiltonians (|30| ) it is sufficient thus to find the matrices 
Mj satisfying (|32D . Using the equations (|26D it is easy to verify that (|32D is satisfied 
with the following matrices |I^: 



Mj{u,X)= I ^_,._^ _^ 1. (33) 

The two properties: canonicity and invariance of Hamiltonians constitute the defi- 



nition of what is called an integrahle map ||3^. It can be considered as a discrete-time 



analog of integrable hamiltonian fiow. Veselov [^ has proved a discrete-time analog of 
Liouville theorem which claims that in the action-angle variables ($, (p) any integrable 
map acts as a shift ipj i— »• ipj + fij($), or, speaking more precisely, as a collection of 
shifts due to the multivaluedness of algebraic mappings. Applying the theorem to the 
case of Backlund transformation depending on a parameter A and noting that shifts on 
the Liouville torus commute we obtain as an immediate consequence the commutativity 
of BT 

Bx, o Bx, = Bx, o Bx,. (34) 

Note that a direct proof of commutativity of BT is not a simple task, see for example 
^ , 15 1 . It is trivialized in our case entirely due to the fact that from the very beginning 



we are working in the hamiltonian context. 

The last property in our list is spectrality. It was discovered rather recently |]19 



and the main motivation for its existence comes from the quantum case, see section 



3.3. 



Let us introduce the quantity fi which is, in a sense, canonically conjugated to A: 

dFx " 

The spectrality of BT means that the pair (e^. A) lies on the spectral curve of the 
Lax matrix. Since det L{u) = 1 it means that both e^ and e~'^ are eigenvalues of L{X) 

Wie^^", A) = det(e±'^ - L(A)) = (36) 

11 



(it does not matter if we take L{X; X, x) or L(A; Y, y) since they are isospectral). 

The property of spectrahty of BT still remains somewhat mysterious and certainly 
needs more research to uncover its algebraic and geometric meaning. The main draw- 
back of the present definition is its being formulated in quite noninvariant terms of 
generating function Fx- 

To prove (^) it suffices to show that, say e^ is an eigenvalue of the matrix L(A; Y, y). 
We shall construct explicitely the corresponding eigenvector uji. 

L{X;Y,y)uJ, = e^uJ,. (37) 

From (|33D it follows that det{Mj{u, X)) = \ — u. It is easy to see that for u = A the 
matrix Mj{X, A) degenerates into a projector 

M,(A,A)=(^'^'^) ^'~'"' ~^^ (38) 

and, as a consequence, has the unique, up to a scalar factor, null-vector 



Un 



Mj{X, X)ujj = 0. (39) 



V W 
Using the identity ( |3TD with M = Mi we conclude that 

Mi{X,X)L{X;Y,y)uJi = (40) 

which, combined with the uniqueness of the null- vector uJi of Mi, implies that ooi is 
an eigenvector of L{X; Y,y). To determine the corresponding eigenvalue, we apply the 
same argument to the identity ( |5^ obtaining the equality Mj+i(A, A)^j(A; Yj, yj)ujj = 
from which it follows that ij{X; Yj,yj)ujj ~ ujj+i. The direct calculation shows that 

i,iX;Y,,y,)uj,=e^^-^-y^uj,^i. (41) 

It remains only to use the formulae (^ and (^5|) to arrive finally at (0). 

An alternative variant of the proof, more close to what we shall use in the quantum 
case (see section ^.51 ) is to introduce a gauge transformation with a triangular matrix 

£,=iVr;,£,(A;r„y,)iV,= (^ e"^, e^-- ^ J ' ^^=(0 1 j ^^2) 

(note that cjj coincides with the second column of A^j). 
The result, as expected, is 

t(A) =trC(A)...£i(A) =trC...£i = e'^ + e-^. (43) 

We conclude this section with a remark on using BT for generating solitons which 



is the main application of BT to the integrable nonlinear evolution equations |12, 15 
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We are following here the argument by Gaudin |T3|. Let us apply the Backlund 
transformation (p6D to the vacuum state Xj = Xj = 0. The equations (|^) turn into 

= e-ys +ey'+' -X, (44a) 

Yj=e-y^ +ey' -X. (44b) 

Concentrating on the first equation (the second equation describes the time evolu- 
tion Yj = di/j/dt with respect to the hamiltonian H) we introduce the parametrization: 
A = 2 cosh K, e^" = cosh(Q; + k)/ cosha. The general solution can be now written as 



^y^ ^ cosh(a + k{j + 1)) 
cosh (a + nj) 



(45) 



In case of the infinite lattice, when j G Z, the formula ( |45| ) describes a soliton 
solution. Note, however, that the solution ( ^51) has different asymptotics e^^ -^ e^*^ as 
j -^ ±oo, satisfying thus the boundary conditions different from those for the vacuum 
state. As a result, the energy and values of other integrals of motion for the soliton 
solution differ from those for the vacuum. 

The situation is quite different in the periodic case. The periodicity condition 
i/n+i = Hi can be satisfyed in two ways. The first one leads to the quantization of the 
parameter: ku G ttIZ and is inacceptable if we want to keep A free. Besides, in this 
way we get a complex solution for e^^. Another option is to fix the free parameter a 
by setting a = ±oo, which gives us another vacuum state e^-' = e^'^ having the same 
values of Hamiltonians as the vacuum. 

The fact that BT in the periodic case does not produce solitons and always preserves 
the integrals of motion may dissapoint those accustomed to other usages of BTs. A 
merit of our variant is, however, that it has deep analogies in the quantum case, as we 
shall see further. 

2.5 Duality 

Besides the 2x2 Lax matrix L{u) which we used until now there exists another, n x n 
Lax matrix C{v) for the Toda lattice which is dual to L{u) in the sense that the cor- 
responding spectral curves are equivalent up to interchanging the spectral parameters 
u and V 

(-1)"-^ det(n - Civ)) = det(t; - L{u)). (46) 

Referring the reader to the paper |I| where the geometric meaning of the duality is 
elucidated, we present here a more elementary approach. 

To produce the dual Lax matrix C{v) we take an eigenvector 6i{u) of L{u) corre- 
sponding to the eigenvalue v (for brevity, we will not mark the dependence on u in 

0) 

L{u)ei = v9i (47) 

and define by induction 6j as 

e,+, = i^{u)e„ j = l,...,n. (48) 
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From ( P^ ) it follows that ft„+i = v9i. The function 6j{u), when properly normalized, 
is called Baker- Akhiezer function. Denoting the components of the vector 6j as ipj and 
ipj we write down (|48|) explicitly as 






u + Xj — e^J' 
e-^J 






Then, splitting the components and taking into account the quasiperiodicity con- 
dition 9n+i = v9i we arrive at the following linear equations for ipj and ipj-. 






(49a) 
(49b) 



V'i+i = e ^'^j, j = 1, ...,n-l 

vi/ji = e~''''ipn. 
Eliminating ipj we obtain a second-order finite-difference equation for cpj 



uipi = ip2-Xnpi + e''^ """V ^ipn, 



ULpj = yjj+i — Xj(fj + e' 



^i ^ J - 1 



^j-1, 



j = 2,...,n-l 



V^n-l, 



(50a) 
(50b) 

(51a) 
(51b) 
(51c) 



which can be rewritten as the linear problem for the vector $ with the components cpj 
in the matrix form: 



Civ] 





£{v)<^ -- 


= «$, 




$= ... 


C{v) defined as 








( -X, 


1 

-X2 .. 








y-ie^i„ \ 






.. 
.. 


. -Xn^ 


-1 
1 


1 



•^jk — -^j '^k 



(52) 



is the dual Lax matrix we were looking for. 

We leave the proof of the identity ( |46| ) as an exercise to the reader. For the r-matrix 
corresponding to the Lax matrix C{v) see pO|, p!6| . 

Similarly to the case of 2 x 2 matrix L{u), for C{v) there must also exist a Darboux 
matrix A4 intertwining C{v]X,x) and C{v]Y,y). The explicit expression for Ai, like 
the one for C{v), can be found from the Baker- Akhiezer function. Let dj and 6j refer. 
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respectively, to £(f;X, x) and C{v;Y,y). Let us assume that 9j and 9j are linked by 
the relation 9j = Mj9j, which is obviously compatible with (^51) and (^2]). Expanding 

^j = Mj9j as 

/(/?j\_/m-A + e^^-^^-i -e% \ / ^j \ 

taking its first line 

ifj = {u - \ + e^^"^-'-^)(^j - e^^'V'j 

and substituting u(pj = (fj+i — Yj'^j + e^J~^^"^(^j_i ipj = e~^^^^(^j_i from 6'j+i = 
£j(ti; Yj, yj)9j, as well as Yj = e^^"^^ + e^^"^'^"^ — A from (]26|b), we obtain, after making 
the necessary correction for j = n the following result: 



iPj = iPj+i -&"' y^ifj, 



iPn = Vipi-e''" y-ifin 



j = l,...,n- 1 



(53a) 
(53b) 



or, in matrix form, = A^G, with 



M{v) = 


{ -e^i-s^i 





1 ... 

_QX2-y2 _ _ _ g 

g QXn-l-yn-l 

... 


\ 


1 

ci^n 1/71 

^ 1 


By construction, we have 






M{v)C{v 


Y,y) = Civ;X,x)Miv). 




Alternatively, one could introduce Mj ~ — M~^ 





(54) 



(55) 



Mj(m,A) 



1 -e% 



such that 



Mj+i(u, A)£j(m; Xj, a;^) = ij{u\ Yj, yj)Mj{u, A) 



and repeat the same calculation, starting from 9j = Mj9j. The result is G = j\4Q, 
with 



M{v) 



( 1 




. 





— V' 


-i^yi-xn \ 


-e^2- 


-XI 


1 . 















. 


1 







V 




. 


Qiln—Xn-l 




1 / 



(56) 



satisfying 

Miv)C{v; X, x) = C{v- r, y)M{v) (57) 

Despite the fact that M. ^ M.^^ the formulas (|^) and ( |57D are compatible because 
of the the remarkable factorization of C{v): 

C{v-X,x)-Xl = M{v)M{v), C{v]Y,y)~\l = M{v)M{v), (58) 
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see [^, ^ for discussion of the factorization as a mechanism for generating Backlund 
transformations. 

In the above formulas v is, by definition, an eigenvalue of L{u), so the pair (w, u) lies 
on the spectral curve det(f — L{u)) = of L{u). When dealing with L{u), it is conve- 
nient to take u as independent variable, and when dealing with C{v) respectively v. For 
Backlund transformation it means in fact swapping the roles of A and /i: the parameter 
H becomes independent numeric variable instead of A. All the formulas defining BT re- 
main the same but their interpretation changes: the equality (^) becomes a constraint 
for X and y rather than definition of /i, whereas A becomes a dynamical variable — a 
Lagrange multiplier for the constraint which can be determined from equations (^). 
The respective dual Backlund transformation i3^ possesses all characteristic properties 
of BT which can be proven using the Lax matrix C{v) in the same manner as for Bx, 
see [|I9| for details. 

2.6 General construction of Backlund transformation 

As shown in section p.4|, to any Backlund transformation Bx there corresponds a Dar- 



boux matrix M(u,X) intertwining the corresponding Lax matrices, see formula (pl|). 
In practice, however, one usually does not know the BT apriori, and has to deal with 
the inverse problem: given L{u) how to find admissible M{u, A) producing a BT. If 
one is not interested in the Hamiltonian properties of the transformation the usual 
strategy is to try some ansatz for M{u, A), say, as a low-degree polynomial in u. See 



the monograph |22[ for a plentitude of examples. 

In this section we shall restrict our attention to the integrable models generated 
by the quadratic Poissonian algebra (JT^) with the S'L(2)-invariant r-matrix (^) and 
address the following question: which M{u, A) are admissible that is produce canonical 
mappping i^A? 

Answer: It is sufficient that M{u, A) as a smooth manifold coincide with a sym- 
plectic leave of the same quadratic Poisson bracket (|14D as L{u), the leading coefficients 



M^™) of M{u, A) and L^") of L{u) in u commute: 

[M(™),L(")] = (59) 

and, also M^^^L^""^ ^ (nondegeneracy condition). 
Open problem: Are these conditions necessary? 
The Backlund transformation Bx constructed for the Toda lattice in section |2^ also 



fits our scheme. Indeed, the Darboux matrix M{u, A) given by ( |5BD has, as a smooth 
manifold, the same structure as the local Lax operator (p^ for the DST model. The 
parameter A is introduced through the shift u ^^ u — X which is an automorphism of 
the Poisson algebra. To elaborate, let M{u, A) be 

M(u, A; 5, s) = i°'^{u -X;S,s)=(''~ ^^+ ''^ ^ J J . (60) 



As we shall see, the equation (|3TD allows then to determine S, s and, eventually, 
L{u;Y,y) in terms of L{u]X,x). Expand first (|3TD in powers of u using ([T9| ) and 
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The coefficient at m""*"^ vanishes because of ([55|). The matrix element 21 of the 
coefficient at m" gives the expression for S: 

S = L^,-'\x,x). (61) 

To determine s take again (0) and substitute u = \. Multiplying the resulting 
matrix equality by the row-vector (1, — s) and noting that at -u = A matrix M('u, A) 
degenerates: 

M(A,A)=(^ j) (^ -^) , (62) 

we obtain the quadratic equation for s: 

Li2(A; X, x) + s(Ln(A; X, x) - L22(A; X, x)) - s^L^iX; X, x) = 0. (63) 

Expressing the variables S and s in terms of X and x one can, in principle calculate 
the Poisson brackets for L{u; Y, y) directly and verify that they have the same r-matrix 
form (|l^) as for L{u; X, x) proving thus the canonicity of the transformation from (X, x) 
to {Y, y). See ||3^ where it is done in a slightly more general situation. The calculation 
by brute force, however, is not particularly instructive, and below, following |^, we 
present a quite simple and general proof. The construction we describe is mimics the 
construction of the quantum Q-operator described in section (|3.4|). 



Suppose that M{u) is a symplectic leaf of the same Poisson algebra (14) as L{u): 



{M{u),M{v)} = [ri2{u-v),M{u)M{v)] (64) 



satisfying the condition (^9]). Note that M{u) is by no means restricted to (P^'^{u — 
X;S,s) as above. Let M{u) be parametrized by the canonical variables {S,s), and 
L{u), respectively, by (X, x). The matrix M{u) might contain one or more parameters A 
which we neglect. Assuming the commutativity (|59|) , consider two products: M{u)L(u) 
and L{u)M{u). By virtue of the comultiplication property of the bracket (|1^) they 
both are symplectic leaves of the same bracket. Furthermore, due to the condition 
j^{m)j^{n) _ ]^in)j\^(m) _^ g^ they share the same values of casimirs described in section 
p.3| , namely, the leading coefficient and determinant. In a generic situation, provided 
there is no accidental degeneration, which we shall assume, the equality of casimirs 
implies an isomorphism of the symplectic leaves, or, in other words, there should exist 
a canonical transformation £ : (X, x; S, s) — > [Y, y, T, t), determined from the equation 

M{u; T, t)L{u; Y, y) = L{u; X, x)M{u; S, s). (65) 

Suppose that the canonical transformation £, has a generating function F{t, y; s, x) 

X.?L, y^.f, s^f, r^-f (66) 

OX oy OS at 

(for simplicity, we omit the indices j in Xj, Xj etc.) 
Let us impose now the constraint 

t = s, T = S (67) 
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and note that on the constraint surface we have M{u; T, t) = M{u; S, s), and therefore 
the equahty ([65| ) is transformed to the Darboux form (0). It remains to prove that the 
transformation £ remains canonical after being restricted on the constraint surface. 

Suppose that one can resolve the equations dF/ds + dF/dt = with respect to 
s = t and express X and Y from (^) as functions of {x, y). 

Proposition. The resulting transformation B : {X, x) -^ {Y, y) is canonical and is 
given by the generating function $(x,y) = F{s{x,y),y] s{x,y),x), such that 



X 



9$ 

dx ' 



Y 



dy' 



Proof. Let \st mean the restriction on the constraint manifold s 
proof consists of two lines: 



(68) 
s{x,y). The 



X 




ds dF 
dx dt 



St 



dt_ 

dx 



St 



We observe now that 



'dF dF' 
ds dt 







(69) 



(70) 



St 



due to S = T, and, consequently, X = d^/dx. Similarly, one establishes Y = —d^/dy 
completing thus the proof. 

In many applications the Lax matrix L{u), like in Toda case, is a monodromy 



matrix factorized into the product of local Lax matrices 
the same Poisson brackets ( p!4D as L{u) 



(u), see formula (H), having 



{i,{u)},{v)} = [ru{u-v)},{ufe,{v)]6i,. 



(71) 

The similarity transformation ( |3l]) is replaced now with a gauge transformation 
(|32| ) which ensures the preservation of the spectral invariants of L{u). 

The modification of the reduction procedure described above is quite straightfor- 
ward. Supposing that ij{u) and Mj{u) depend on local canonical variables we define 
first the local canonical transformations £*•*-' : {Xj^xf Sj, Sj) -^ (Yj,yj,Tj,tj) from the 
equations 

Mj{u] Tj, tj)ij{u] Yj, yj) = ij{u; Xj, Xj)Mj{u] Sj, Sj). (72) 

Let the corresponding generating functions be f^^\tj,yj;Sj,Xj). Consider the di- 
rect product of ra phase spaces {Xj,Xj] Sj,Sj) and, respectively, (Yj,yj]Tj,tj). The 
generating function 

n 

F:=Y.f^'\t„yr,s„x,) (73) 

3=1 

determines then the direct product £ of the local canonical transformations £,^-^'. 
Let us now impose the constraint 



U 



Sj+l, 



S, 



(74) 



assuming periodicity j + n = j. The proof of the canonicity of the resulting transfor- 
mation B : {X, x) -^ (Y, y) parallels the proof given previously. It remains to notice 
that after imposing the constraint ([7^) we have Mj{u; Tj, tj) = Mj{u; Sj, Sj) and obtain 
the equality (p2D. 

It is convenient to represent the structure of the Backlund transformation graphi- 
cally. Let the local transformation £,^^' be depicted as a four-legged vertex (see figure 
|l]), each leg corresponding to a canonical pair like {X,x) etc. The arrows show the 
direction of the transformation. 



Vj 






Figure 1: Local transformation 

The Backlund transformation B is represented then by the figure @ where the joint 
horizontal lines mark the constraints (73). 



^n t-n— 1 

< 



Vn-l 



1/2 



S2 = h 

* 



yi 



Si tfi 

< ■ 



Xn 3^71—1 -^2 -^1 

Figure 2: Composition of local transformations 

In conclusion to this section, a few general remarks. The proof of canonicity pre- 
sented above is pretty general using only the comultiphcation property of the quadratic 
r-matrix Poisson bracket ([T^). It covers thus all integrable models governed by the 
bracket ( jT^ ) for any r-matrix, not necessarily the 5'L(2)-invariant one. The only thing 
one needs is to study the structure of symplectic leaves of the Poisson bracket and to 
choose some elementary matrices M{u). 

Note that the product of two M- matrices produces the composition of corresponding 
Backlund transformations. Given the conjecture about the factorization of symplectic 



leaves from section p.3| is true, it implies that any BT is decomposable into elementary 



BTs corresponding to above the Lax matrices 



[u 



\),i- 



[U 



X),i- 



[U 



A), 



r°''^(M-A),r°'^^(M-A). 

An interesting and yet unsolved problem is how to deal with the spectrality property 
of BT within our construction. Our conjecture is that there is a spectrality identity 
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det(e'^ — -^(A)) = with respect to any zero m = A of det M{u). 

2.7 Application to Toda lattice 

Let us demonstrate how the construction described in the previous section produces 
the Backlund transformation for the Toda lattice described in section |2.4| . Substituting 
into the formula (|72D the expressions ( p!0D for ij{u) and (|60D for Mj{u, A) we get 



J 



f U-X + tjTj -tj \ ( u + Yj -ey^ \ _ ( u + Xj -e^^ \ ( u-X + SjSj -s 
[ T, -1 ) [ e-y^ )-[ e'-^ )[ S, -I ) 

(75) 
The system of equations obtained by equating the coefficients at powers of u has a 
unique solution: 

Yj = -A + e^^sj^ + SjSj, (76a) 

ey^=Sj, (76b) 

Tj=e-''\ (76c) 

tj = Xe^^ - e^'^^sj^ + e^^Xj. (76d) 

defining the local transformation £^ . Strictly speaking, due to the degeneracy of 
the Lax matrix i'^°'^'^{u) the proof of the canonicity of Sl^ given in section ( p.6| ) does 
not apply here directly because the transformation ([76| ) does not possess a generat- 
ing function in terms of (t,y,x,s). It is easy, however, to verify the canonicity by a 
straitforward calculation. 

The equalities (|76|b) and (|76|c) allow to resolve the constraint (0) yielding Sj = e^^' , 



Sj = e~^^'^ which, upon being substituted into (|76| a) and (|76|d) produce exactly the 
defining relations ( p6D for the Backlund transformation studied in the section p.4| . 

Exercise. Find what canonical transformations preserving the Hamiltonians of the 
Toda lattice are generated by the Darboux matrices M{u) = diag(l, a) and M{u) = 

r°'^^(u-A). 

3 Quantization 

3.1 Quantum/classical correspondence 

Here we give only a very brief account of the quantum mechanical notions we are 
going to use. For more information on the basics of quantum mechanics see any good 
textbook. See also Reshetikhin's lectures in this volume for references on deformation 
quantization. 

The quantum observables are usually introduced as self-adjoint operators in a 
Hilbert space. In the limit of the classical Hamiltonian mechanics, as the deforma- 
tion parameter h (Planck constant) goes to 0, the observables commute, and the next 
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order term in h produces the Poisson bracket of the corresponding classical observables: 

[.,.] = -ih{.,.} + Oih'). 

We shall work with the realization of the Hilbert space of quantum states as a 
space £2(K") of square integrable functions of canonical coordinates (xi, . . . , x„). The 
corresponding canonical momenta are quantized then as the differentiation operators 
Xj = —ihd/dxj. Generally speaking, any operator Q in L2{M."') can be realized as an 
integral operator 

Q : fix) ^ j dxi... dxn Q{y I x)f{x) (77) 

with the kernel Q{y \ x) which possibly is a generalized function (distribution). 

To the canonical transformations in the classical mechanics (automorphisms of the 
Poisson algebra) in the quantum mechanics there correspond the automorphisms of 
the associative operator algebra, that is similarity transformations A i— >■ QAQ~^ with 
unitary operators Q. The following beautiful formula [^ 

Q(y I x) ~exp(i^-^F(y I x)), /i -^ (78) 

gives the correspondence between the kernel Q{y \ x) oi a. unitary transformation and 
the generating function of the classical canonical transformation into which it turns in 
the classical limit. The formula (|78|) works for non-unitary transformations as well. 

In what follows we shall occasionally use non-self-adjoint and non-unitary opera- 
tors which corresponds in the classical case to working with complex rather than real 
manifolds. 

3.2 Quantum Toda lattice 

Starting from now we shall drop h from our formulas assuming h = 1. The commutative 
Hamiltonians of the periodic quantum Toda lattice are differential operators in £2(1^")- 
They are obtained from exactly the same formulas (|]), (P), ([10|), ([l^) as in the classical 



case where one should substitute Xj = —idx- The proof of their commutativity is 
based on the algebraic framework called The Quantum Inverse Scattering Method, see 



18, 31, 32] for a detailed exposition of the method. 



Starting with the quantum local Lax matrix £,• (u) 




we observe that it satisfies a quadratic commutation relation 

Ruiui - U2yiuifiiu2) = }iu2fiiui)Ri2iui - U2) (79) 

with 

Ruiu) = u + iVi2. (80) 
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As its classical counterpart (JT^), the relation ( |75D possesses the comultiplication 
property which implies that the monodromy matrix L{u) given by (|^) satifies the same 
relation 

12 2 1 

Ruiui - U2)L{ui)L{u2) = L{u2)L{ui)Ri2{ui - U2) (81) 

from which the commutativity of the Hamiltonians 

[t(Mi),t(M2)]=0 



follows immediately (see [|1^, |3T|, |32| for explanations). 



The associative algebra given by the generators L{u) and quadratic relations ( pT]) 
is called yangian 3^[5'/2]- Its representations correspond in the classical limit to the 
symplectic leaves of the quadratic Poisson bracket (jT^). A convenient way of viewing 
the equality (0) is to treat it as a particular form of the quantum Yang-Baxter equation 

/2i2(m)/?i3(m + V)R2M = R23iv)Rl3iu + v)R2siv) (82) 

which is considered as an operator equality in the tensor product V^i (8> V2 ® V3 of three 

linear spaces Vi, V2 and V3. Respectively, Rjkiu) is an operator in the space Vj ®Vk 

naturally embedded in Vi®V2®V^. For any pair Vj, Vk of the yangian 3^[5'/2] nioduli 

there exists Rjk such that the Yang-Baxter equation (|82D holds for any triplet V^i, V2, 

V3. 

In particular, for Vi = V2 = V3 = C^ we have the YBE (|^) for the i?-matrix (0). 

The relation (^) also can be considered as a particular case of (|82D for Vi = V2 = C^ 

1 2 

(auxiliary spaces), V3 = £2(1^") (quantum space), i?i3 = L, R2Z = L. 

The commutativity of the quantum Hamiltonians being established, the next prob- 
lem is to find an effective way of determining their joint spectrum. There are two known 



ways of approaching this problem: separation of variables pl| , p6| , |17| and Q-operator 
I^ p5| . Here we shall consider the latter approach. 



3.3 Properties of Q-operator 

The original idea of Baxter [|^, |^ which enabled him to solve the XYZ spin chain was to 
construct a one-parametric family of operators Q\ commuting with the Hamiltonians 
of the model 

[QA,t(n)]=0 (83) 

and hence sharing with t{u) the common set of eigenvectors. Moreover, Q\ must satisfy 
the Baxter equation 

QxtiX) = A+(A)g,^i + A_(A)g,_i (84) 

where A±(A) are scalar functions determined by the parameters of model. Note that 
in the left-hand-side of ( ^4]) the order Qt or tQ is not important because of the com- 
mutativity (^3|) . Applying the Baxter equation (^) to a common eigenvector of Q\ 
and t(A) one can replace the operators in (^^ by their eigenvalues. The resulting 
finite-difference equation of second order for the eigenvalues of Qx considered in an 
appropriate functional class allows then to determine the spectrum of t{X). 
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Baxter succeeded to construct a Q-operator for the XYZ spin chain as a trace of a 
monodromy matrix 

gA = tryL„(A)...Li(A) (85) 

constructed with a specially chosen auxiliary space V. Graphically, the structure of Qx 
is represented by the samefigure |^ which we used in the classical case. The horizontal 
lines correspond the auxiliary space, the vertical ones — to the quantum space. Each 
vertex represents an operator hj{X). The commutativity (|83|) is garanteed then by 
the Yang-Baxter equation, and the only problem is to choose such V which would 
produce the Baxter equation (|8^). Later Gaudin and Pasquier [^ have constructed a 
Q-operator for the quantum periodic Toda lattice by giving an explicit expression for 
its kernel Qx{y \ x) as an integral operator. They have also noticed that the classical 
limit of the similarity transformation Qx{-)Qx^ is exactly the Backlund transformation 
studied in the previous sections. 

Below we reproduce the result by Gaudin and Pasquier. Our approach combines 
their integral operators technique with Baxter's original idea of constructing as a trace 
of a monodromy matrix (|85D . 

Note that such properties of Backlund transformation as the invariance of Hamil- 
tonians (|30|) and spectrality (|43|) are the classical counterparts of such properties of 



Q-operator as, respectively, commutativity (|^) and the Baxter equation (p^ . The 
former one being obvious, we comment only the latter one. Observing that the shift 
operators A t— > A =F i are expressed as exp(=i=i9A) = exp(±/i) where /i is the canonical 
momentum conjugate to A we can rewrite (Q) in the form t(A) = A_|_(A)e~'^ + A_(A)e^ 



which gives (|43| ) in the classical limit (for the Toda lattice Aj. = 1) 



3.4 Q-operator for Toda lattice 

We shall construct the Q-operator as the trace of the monodromy matrix (^) taking 
for the auxiliary space V the space C[s] of polynomials in variable s. The corresponding 
representation of the yangian 3^ [(7/2] is realized then as the Lax operator of the quantum 
DST model 

7. ^/ A \ I u — X — ids —s\ /„„N 

M{u,X)= _i5 _i , (86) 



compare with the formula ( |5DD for the classical case. 

To prove the commutativity (|83| ) it is sufficient to establish the identity 



M{u, X)i{u)hx = hx£iu)M{u, A) (87) 

which can be considered as a variant of the YBE (|8^) with the following layout of 
spaces: Vi = C^ V2 = C[s], V3 = £2(^")- We shall use (0) as the equation for 
determining L^. Rewriting (^) as the system of equations for the kernel hxit, y \ s,x) 
of La 

u — X — itdt —t \ I u — idy — e^ \ ^ /, 1 n 

'J Q j^x{t,y\s,x) 

.^^ _^ j^x{t,y\ s,x) (88) 
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we obtain a unique, up to a scalar factor, solution 

hx{t, y\s,x)^ 6{s - e^) exp(ite-^ - ie^"^ + iA(x - y)). (89) 

From (^) we get the formula for the kernel of Qx: 

Qx{y \ x) = dsn... dsi W La(sj+i, yj \ Sj,Xj). (90) 

i=i 

The integration over Sj in (pO| ) reduces, due to the delta-function factor in (|89|), to 
the substitution Sj = e^^. Finally, we have 

n 

Q^{y I x) = n exp(ie^^+i-^^ - ie^^"^^ + iX{xj - yj)). (91) 

j=i 

Note that Qx{y \ x) = exp {—iFx{y \ x)) where Fx{y \ x) is the generating function 
( |29| ) of the classical BT, that is the semiclassical formula (^) is exact in our case. 
This is an accidental peculiarity of Toda lattice which usually does not hold for other 
models. 



In [^ another version of the Q-operator is used which differs from (|9lD by the shift 



Vj ^ Vj + i^/2 



Qx{y \x) = Y[ exp(-e^^+i-^^ - e^^-y^ + X{x, - y^)) (92) 



which, in operator terms, corresponds to multiplying Qx by the factor exp(— 7ri7i/2). 
The kernel (|92D is more convenient for analytical study since it rapidly decreases along 



the real axis in Xj. 



3.5 Baxter's equation 

The commutativity (|83D being already established, it remains to prove for our Qx the 
Baxter equation (0). We reproduce here the proof by Gaudin and Pasquier [|25| which 



parallels the proof for the classical case given in the end of section p.4| , see formula 

(El)- 

First, note that the kernel (0) factorizes as 

n 

Qxiy\x) = l[^jW (93) 

into factors 

Wj{X) = expiiey^-""^-' - ie^'^-y^ + iA(a;j_i - yj)) (94) 

Applying then t(A) to the kernel Qx{y \ x) and using (^ we observe that each 
ij{\; —idy.,yj) acts locally only on Wj{\) and obtain 

t{\)Qx{y I x) = tr(4(A)w;„(A)) . . . {h{\)wi{\)) = Qx{y | x) tr£„ . . . £i (95) 
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where 

£,=£,(A)lni^,(A)=(^" /^; I y (96) 

After that we can use the triangular gauge transformation ij = N'^iijNj with Nj 
and the resulting matrix ij given by the same formulas (^2]) as in the classical case. 
Noticing then that 

w,{X) ' w,{X) ^"^ 

we obtain the required result 

t{\)Qx = Q^^i + Q^_i 

Similarly, for the modified kernel ( ]92| ) one obtains 

t{X)Qx = i"Q,+i + i-"Q,„i. (98) 

The Toda Hamiltonians {Hj}"'^^ enter the Baxter equation (|98D through the gen- 
erating function t{u) = m" + Hiu"'~^ + . . . + if„. Their eigenvalues are determined by 



the condition that the finite-difference equation (|98D possesses a solution Qx which is 
holomorphic and rapidly decreases along the real axis. For a detailed analysis of the 
equation (|98D see ||2^, [1] 



4 Conclusion 

We have discussed here only the most elementary properties of Backlund transforma- 
tion and Q-operator using the sole example of Toda lattice. For further reading see 
], 0,|21|, |2|, H, 13. 
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